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ABSTRACT 

The origin of ultra-intense magnetic fields on magnetars is a mystery in modern astro- 
physics. We model the core collapse dynamics of massive progenitor stars with high 
surface magnetic fields in the theoretical framework of a self-similar general polytropic 
magnetofluid under the self-gravity with a quasi-spherical symmetry. With the speci- 
fication of physical parameters such as mass density, temperature, magnetic field and 
wind mass loss rate on the progenitor stellar surface and the considcration of a re- 
bound shock breaking through the stellar interior and envelope, we find a remnant 
compact object (i.e. neutron star) left behind at the centre with a radius of ~ fO 6 
cm and a mass range of ~ f — 3 M Q . Moreover, we find that surface magnetic fields 
of such kind of compact objects can be ~ 10 14 — fO 15 G, consistent with those in- 
ferred for magnetars which include soft gamma-ray repeaters (SGRs) and anomalous 
X-ray pulsars (AXPs). The magnetic field enhancement factor critically depends on 
the self-similar scaling index n, which also dctcrmincs thc initial density distribution 
of the massive progenitor. We propose that magnetized massive stars as magnetar 
progenitors based on the magnetohydrodynamic evolution of the gravitational core 
collapse and rebound shock. Our physical mechanism, which does not necessarily re- 
quire ad hoc dynamo amplification within a fast spinning neutron star, favours the 
'fossil field' scenario of forming magnetars from the strongly magnetized core collapse 
inside massive progenitor stars. With a range of surface magnetic field strengths over 
massive progenitor stars, our scenario allows a continuum of magnetic field strengths 
from pulsars to magnetars. The intense Lorentz force inside a magnetar may break 
the crust of a neutron star into pieces to various extents. Coupled with the magne- 
tar spin, the magnctosphcric configuration of a magnetar is most likely variable in 
the presence of exposed convection, differential rotation, equatorial bulge, bursts of 
interior magnetic flux ropes as well as rearrangemcnt of broken pieces of the crust. 
Sporadic and violent releases of accumulated magnetic energies and a broken crust are 
the underlying causes for various observed high-energy activities of magnetars. 

Key words: magnetohydrodynamics (MHD) — shock waves — stars: magnetic fields 
— stars: neutron — supernova remnants — white dwarfs 



1 INTRODUCTION 

Magnetars are believed to be neutron stars with surface 
magnetic field strengths considerably stronger than the 
quantum critical value of Bqed = 4.4 x 10 13 G. There are 
two main types of observational manifestations for magne- 
tars: (i) Soft Gamma-ray Repeaters (SGRs) and (ii) Anoma- 
lous X-ray Pulsars (AXPs). Up to now, six SGRs and ten 
AXPs have been identified observationally (see Mereghetti 
2008 for a latest list and an extensive review as well as very 
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recent powerful explosions of SGR J1550-5418 with a short- 
est spin period of 2.07 s). Most recently, a new Galactic mag- 
netar is reported with very fast optical flares (|Kouveliotoul 
l2008l ; IStefanescu et al.ll2008l ; ICastro-Tirado et al.ll2008l ). al- 
luding a continuum from ordinary dim isolated neutron stars 
to magnetars. The ultra-intense surface magnetic fields on 
magnetars are unique in the Universe and they are respon- 
sible for various high-energy activities, for example the gi- 
ant 7-ray fla re of SGR 1806-20 (e.g. iHurlev et al] 120051 ; 
iPalmer et ai1l2005h . Magnetar-like X-ray emissions are also 
detected from a rotation-powered pulsar PSR J1846-0258 
with an inferred intense mag netic field of ~ 4.9 x 10 13 G at 
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the centre of supernova rem nant Kes75 (e.g. iGavriil et al.l 
l200Sl ; lArchibald et aLlliooih . 

Recent observations have also provided clues connecting 
magnetars with very massive progenitor stars, for example 
an infrared elliptical ring or shell was discovered surround - 
ing the magnetar SGR 1900+14 (e.g. IWachter et aljlioosh . 
However, the formation of magnetars, especially the origin of 
the ultra-intense magnetic field, remains an important open 
issue. There are two major contending physical scenarios, 
viz. the dynamo scenario versus the f ossil-field scenario. 

Duncan fc Thompsonl (|l992h and lThompson fc Duncanl 



couple of early B stars, for example th e B0.5V star HD 370 61 
(~ 650 G, e.g. Hubrig et al. 2006). iPetit et al.l (|2008al lbh 
carried out systematic spectropolarimetric observations to 
search for magnetic fields on all massive OB stars in the 
Orion Nebula Cluster star-forming region. Strong magnetic 
fields of the order of kG were inferred on 3 stars out of a 
sample of 8. The existence of strong magnetic fields on OB 
stars even appears somewhat overwhelming in contrast to 
very few magnetars that have been discovered so far. 

With the assumption that neutron stars form dur- 
ing the collapse of massive progenitors in the Galac- 
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ring primarily in the convection zone of the progenitor, as 
well as in a differentially rotating nascent neutron star, and 
concluded that very strong magnetic field, in principle up 
to ~ 3 x 10 lr G, may be created. The dynamo mechanism 
requires an extremely rapid rotation of a nascent neutron 
star with a spin period of a few milliseconds. However, the 
current population of magnetars appears to be slow rota- 
tors, having spin periods in the range of ~ 2 — 12 s (e.g. 
lMereghettill2008h . Therefore, neutron star dynamo scenario 
for magnetars faces a considerable challenge to account for 
the fact of slowly rotating magnetars as observed so far. 

The fossil-field scenario for the mag- 
netism of compact objects was first pro- 
posed to explain mag ne tic white dwarfs (e.g. 
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2005; iFerrario fc Wickram asinghc 2005; Lou fc Wang||2007h . 
It is conceivable that the magnetic field of white dwarfs may 
be of fossil origin from the main-sequence phase of their 
progenitors, and the attempt to link magnetic white dwarfs 
with their main-sequence progenitors naturally makes the 
chemically peculiar Ap and Bp stars as pl ausible candi- 
dates. Observations of lAuriere et al.1 l|2003h have shown 
that chemically peculiar Ap and Bp stars are generally 
magnetic indeed, with a surface magnetic field of ~ 100 G 
by Zeeman splittings. In general, magnetic field strengths 
fall in the range of ~ 3 x 10 2 - 3 x 10 4 G (e.g. Braithwaite & 
Spruit 2004 and references therein). Magnetic white dwarfs 
may be created as a result of rebound shock explosion 
l|Lou fc Wand [2007) and may further give rise to novel 
magnetic modes of global stellar oscillations (Lou 1995) . By 
the magnetic flux conservati o n duri ng the stellar evolution, 
iFerrario fc Wickramasinghel (|2005l ) argued that stellar 
magnetic fields (~ 100 G) in their main-sequence phase 
can be enhanced up to the range of ~ 10 6 - 10 9 G on the 
surface of magnetic white dwarfs. This fossil-field scenario 
is supported by the statistics for the mass and magnetic 
field distributions of magnetic white dwarfs. 

B ased on the same scenario an d a sim ilar physical argu- 
ment. iFerrario fc Wickramasinghel (|2006h further suggested 
that the ultra-intense magnetic field over the surface of mag- 
netars may also come from fossil magnetic fields. The pro- 
genitors of magnetars are expected to be O stars and early 
B stars with high surface magne tic fields of ~ 1000 G. 
iFerrario fc Wickramasinghel (|2008h presented a population 
synthesis study of the observed properties of magnetars and 
found that magnetars arise from high mass progenitor stars 
(2OM < M < 45M ). 

Up to now magnetic fields are directly measured for two 
O stars, namely ^Ori C (~ 1 kG, e.g. Donati et al. 2002) 
and HD 191612 (~ 1.5 kG, e.g. Donati et al. 2006), and a 



the main-sequence phase), and ~ 8 percent of massive 
st ars have surface magnetic fi elds higher than ~ 1000 
G, iFerrario fc Wickramas inghc (200(1) estimated that these 
high-ficld massive progenitors gave birth to 24 neutron stars 
with magnetic field > 10 14 G, consisting a major part of 
magnetars. While the fossil-field scenario appears plausible 
from the perspective of statistics, it is highly instructive to 
have a more direct magnetohydrodynamic (MHD) model 
description for the core collapse of high-field massive pro- 
genitor stars and to check whether compact remnants left 
behind MHD rebound shocks do possess ultra-intense mag- 
netic fields. 

In this paper, we attempt to model magnetized massive 
progenitor stars with a quasi-spherical general polytropic 
magnetofluid under the self-gravity (Wang & Lou 2008; Lou 

6 Hu 2009). We examine semi-analytic and numerical so- 
lutions to explore the self-similar MHD evolution emerging 
from dynamic processes of core collapse and rebound shock 
travelling in the stellar envelope with a wind mass loss. More 
specifically, we adopt a general polytropic equation of state 
(EoS) p — n(r, i)p 7 with p, p, 7, and k respectively be- 
ing the gas pressure, mass density, polytropic index and a 
proportional coefficient dependent on radius r and time t. 
Here, k is closely related to the 'specific entropy' and is not 
necessarily a global constant. By 'specific entropy' conser- 
vation along streamlines, another key parame ter q arises in 
self-similar dynamics (see IWang fc Loul l2008h . For being 
a global constant, or equivalently q = 0, the general poly- 
tropic EoS reduces to a conventional polytropic EoS. By 
further setting 7 = 1, a conventional polytropic gas reduces 
to an isothermal gas (e.g. Lou & Shen 2004) . We also require 

7 ^ 1 to ensure a po s itive specific enthalpy p/(~f — 1). 

IChiueh fc Choul l|l994h studied the isothermal MHD by 
including the magnet ic pressure gradi ent force in the radial 
momen tum equation. lYu fc Loul (120051). Yu. Lou, Bia n & Wu 
(2006), IWang fe Loul (f2007h and l Wang fe Loul (|2008h gener- 
alized the self-similar hydrodynamic framework by including 
a completely random transverse magnetic field with the ap- 
proximation of a 'quasi-spherical' symmetry (e.g. Zel'dovich 
& Novikov 1971); the radial component of such magnetic 
field is much weaker than the transverse components. We 
conceive a simple 'ball of thread' scenario for random mag- 
netic fields in a massive progenitor star. In other words, a 
magnetic field line follows the 'thread' meandering within 
a thin spherical 'layer' in space in a completely random 
manner. Strictly speaking, there is always a weak radial 
magnetic field component such that field lines in adjacent 
'layers' remain physically connected throughout in space. 
In particular, we emphasize that the nature of the ran- 
dom magnetic fields inside a progenitor star may be ei- 
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ther those generated by dynamo mechanism pr obably linked 
with convection and differentia l rotation (e.g. ISpruij [20021 ; 
iHeeer. Wooslev fc Spruidl2005l ). or those of 'fossil fields' en- 
trained from molecular clouds in dynam ic processes of star 
formation (e.g. IWang fc Loull2007l . fSooih . 

According to numerical simulations of differentially ro- 
tating magnetized stars by Heger et al. (2005), the dynamo- 
powered radial magnetic fields in a progenitor star are about 
3 to 4 orders of magnitude weaker than the transverse mag- 
netic fields during the pre-supernova evolution. In line with 
this hint, we simplify the treatment by only dealing with the 
dominant transverse field and focus on their dynamic effects 
on the bulk motion of gas in the radial direction. By taking 
the ensemble average of magnetic fields in each thin spherical 
'layer', we smooth out small-scale magnetic field structures 
and are left with 'layers' of large-scale transverse magnetic 
fields. We also presume that small-scale non-spherical flows 
as a result of the magnetic tension force may be neglected 
as compared to large-scale mean radial bulk flow motions. 
Therefore on large scales, a completely random magnetic 
field contributes to the dynamics in the form of the average 
magnetic pressure gradient force and the average magnetic 
tension force in the radial direction. 

This theoretical model framework of self-similar MHD 
has been applied to gravitational core collapse and rebound 
shock proces s es w i thin progenitor stars for supernovae. 
ILou fc Wand (|2006l . 120071 ) modelled the hydrodynamic and 
MHD reboun d shocks of s uperno vae in the self-similar phase 
of evolution. iHu fc Loul (|2008bh presented preliminary re- 
sults for a shock breakout to reprodu ce the early X-ray light 
curve of supernova SN 2008D (e.g. ISoderberg et alj 120081 ; 
iMazzali et aLll2008I ). In this paper, we demonstrate that such 
a self-similar MHD process may give birth to a compact 
remnant with a nuclear density and a range of ultra-intense 
surface magnetic fields. We will see that massive progenitor 
stars whose collapsing cores have magnetic fluxes similar 
to those of magnetars will eventually collapse into neutron 
stars with a magnetar level of magnetic fluxes because of 
the magnetic flux conservation. In our model scenario, the 
neutron star dynamo processes and the required initial rapid 
spins of nascent neutron stars may not be necessary. 



where u is the bulk radial flow speed and < B 2 > is the en- 
semble mean square of a random transverse magnetic field. 
The weak radial component of the magnetic field is deter- 
mined by equations (10) and (II) of Yu & Lou (2005). With 
the self-similar MHD transformation of Wang & Lou (2008) , 
the ideal MHD equations together with the magnetic flux 
frozen-in condition and the general polytropic EoS can be 
readily reduced to a set of nonlinear MHD ODEs in the 
highly compact form of 

X(x, a, v)a = A(x, a, v) , X(x,a,v)v' = V(x,a,v) , (2) 

where the prime ' stands for the first derivative d/dx, and 
the three functionals X, A and V are defined by 



X(x, a, v) = C 
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(3) 



In this straightforward yet somewhat tedious derivation, we 
have useful relations m = ax (nx — v) and q = 2(n + 7 — 
2)/(3n — 2) and we have performed the following MHD self- 
similar transformation for a general polytropic gas, viz. 

r = k 1/2 t n x, u = k 1/2 t n -\ , 



2 GENERAL POLYTROPIC SELF-SIMILAR 
MAGNETOHYDRODYNAMICS 

Under the approximation of quasi-spherical symmetry and 
based on the physical idea outlined in the introduction, 
the ideal MHD equations involve mass conservation, radial 
momentum equation, specific entropy conservation along 
streamlines and the magnetic induction equation (see Wang 
& Lou 2008). We highlight the essential parts of this formu- 
lation of nonlinear MHD equations below. 



2.1 Theoretical MHD Model Formulation 

The ideal magnetic induction equation (without the resistiv- 
ity) implying the frozen-in condition for the magnetic flux 
can be cast into the following form 



d 
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(3n - 2)G ' 

where G = 6.67 x I0 -8 g _1 cm 3 s~ 2 is the gravitational con- 
stant, M is the enclosed mass at time t within radius r, x 
is the independent self-similar variable, v(x) is the reduced 
flow speed, a(x) is the reduced mass density, m(x) is the re- 
duced enclosed mass, k, n and C are three parameters, and 
the dimensionless coefficient h is referred to as the magnetic 
parameter such that < Bf >= lQh r K 2 Gp 2 r 2 . It follows that 
q = 2/3 leads to 7 = 4/3 for a relativistically hot gas. Only 
for the case of q = 2/3, can the parameter C be indepen- 
dently chosen; otherwise for q 7^ 2/3, C can be set to I with- 
out loss of generality (Lou & Cao 2008; Lou & Hu 2009). The 
magnetosonic critical curve (MCC) is determined by the si- 
multaneous vanishing of the numerator and denominator on 
the right-hand sides (RHS) of ODEs (0 and Once so- 
lutions are obtained for v(x) and a(x), the mean magnetic 
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field strength < B 2 > 1/2 can be readily determined from 
self-similar MHD transformation ((4]). With proper asymp- 
totic conditions as well as eigen-derivatives across the MCC, 
nonlinear MHD ODEs <(2j) and ((3} can be numerically in- 
tegrated by using the standa rd fourth-order Runge-Kutta 
scheme (e.g. iPress et al.lll98rj ). 

It is straightforward to treat MHD shock conditions 
for self-similar solutions to cross the magnetosonic singu- 
lar surface X(x, a, v) = 0. By the conservations of mass, 
radial momentum, and energy as well as magnetic induc- 
tion equation in the comoving framework of reference across 
an MHD shock front, we obtain a set of jump conditions 
for an MHD shock that can be cast into a self-similar form 
(see AppendixQg Wang & Lou 2008 and Lou & Hu 2009). 
Note that the so-called 'sound' parameter k in self-similar 
MHD transformation Q is related to the polytropic sound 
speed and changes across a shock front, with the relations 
k-2 = A 2 fci, hi = fi2 = h, x\ — Xx^ where subscripts 1 and 
2 refer to the immediate upstream and downstream sides of 
a shock and A is a dimensionless scaling parameter. Strictly 
speaking, magnetic fields have very weak radial components 
normal to the shock front. Our treatment of magnetic field 
coplanar with the shock front represents a very good ap- 
proximation for our purposes. 



2.2 Analytic Asymptotic MHD Solutions 

The analytic asymptotic solutions at large x of nonlinear 
coupled MHD ODEs © and © are 
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where A and B are two integration constants, referred to as 
the mass and velocity parameters (|Wang fc Loul [20081) . To 
ensure the validity of asymptotic MHD solution ((jjj, we re- 
quire 2/3 < n < 2; the inequality n > 2/3 is directly related 
to self-similar MHD transformation ((4]) where a positive M 
is mandatory on the ground of physics. For 2/3 < n < 1, 
both mass and velocity parameters A > and B are fairly 
arbitrary. In case of 1 < n < 2, we should require B = 
to avoid a divergent v(x) at large x unless a flow system 
is truncated. Using this asymptotic solution at large x as 
initial conditions, the key scaling index n determines the 
initial mass density distribution as p oc r~ 2//n . In our self- 
similar scenario, the valid range of exponent n corresponds 
to a range of density power laws p oc r~ 3 to p oc r~ L . 

By setting v = for all x in nonlinear MHD ODEs 
(0 and ©, we readily obtain an exact global magnetostatic 
solution, namely 

-2/n 



a = Aox 

where the proportional coefncient Aq is given by 
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n 2 - 2(1 - n)(3n - 2)h 



-1 -l/(n-3n ? /2) 



2(2-n)(3n-2) 
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(7) 



This describes a magnetostatic singular polytropic sphere 
(SPS) with a substantial generalization of q 7^ 0; the case of 



g = 0orn + 7 = 2is included here and corresponds to a 
conventional polytropic gas in magnetostatic equilibrium. A 
further special case of n = 7 = 1 corresponds to a magneto- 
static singular isothermal sphere (SIS). Physically, expres- 
sion (0 requires q ^ 2/3 and h < h c = n 2 /[2(l — n)(3n — 2)] 
for the existence of the global magnetostatic SPS solution 
in a general polytropic gas. For n = 4/5, h c reaches the 
minimum value h c = 4. This places a constraint only when 
n < 1; while for n ^ 1, parameter h > is fairly arbitrary. 

There exists an analytic asymptotic MHD solution ap- 
proaching the magnetostatic SPS solution at small x (re- 
ferred to as the 'quasi-magnetostatic' asymptotic solution), 
namely v = Lx K and a = A x~ 2/n + iV:r A '" 1 " 2/n where K 
is the root of the following quadratic equation 



n 2 j (3n - 2) 
+ nh+ -( 



2(3n - 2) 



K 2 + ^z£L K 



2(2-n)(l-n) n 2 + (3n - 2) 2 (1 - 4/n)Q 



-h + 



(3n - 2) 



, (8) 



where Q = g{n 2 /[2(2 - n)(3n - 2)] - (1 - n)h/(2 - n)} is 
introduced for notational clarity (Lou & Wang 2006, 2007). 
In a certain regime, at least one root of quadratic equation 
((8]) satisfies "l(K) > 1 and therefore quasi-magnetostatic 
solutions do exist. The two coefficients L and N are simply 
related by the following algebraic expression 



i(K - 1)N = (K + 2 - 2/n)A L 



(9) 



In this case, the magnetic Lorentz force (i.e. magnetic pres- 
sure and tension forces together) and the gas pressure force 
are in the same order of magnitude in the regime of small x. 

It can be proved that the parameter regimes where 
quasi-magnetostatic solutions exist are 7 ^ 1, h < h c , 
q < 2/3, n < 0.8 and 7 > 1, h < h c , q > 2/3. With pa- 
rameters outside these two regimes, the so-called strong-field 
asymptotic solutions at small x, for which the magnetic force 
dominat es over the thermal pressure force, have b een shown 
to exist (|Yu et al.ll2006l ; IWang fc Loull2007l . l2008t ) . 



3 FORMATION OF COMPACT MAGNETARS 

3.1 Model Progenitors and Compact Remnants 

To be specific, the radial range of our model solutions is 
set within r; < r < r Q , where r\ = 10 6 cm if the compact 
object is a neutron star or black hole, and r = 10 12 cm 
as the radius of a typical massive main-sequence st ar (e.g. 
iHerrero etail 1 19921 ; ISchonberner fc Harmaned Il995l ) . Mas- 
sive stars may have undergone tremendous mass losses be- 
fore the onset of gravitational core collapses, and a progen- 
itor immediately before a core collapse may have already 
lost the entire hydrogen envelope and become a compact 
Wolf-Rayet star with a radius of ~ 10 11 cm and a mass of 
~ 4 — 8Mq (see e.g. Soderberg et al. 2008 and Mazzali et al. 
2008 for recent observations) . We attempt to further identify 
plausible conditions of forming magnetars by gravitational 
core collapses of massive progenitor stars. 

At the initial time of our analysis, the model should 
approximately represent the final stage of a massive pro- 
genitor star before the gravitational core collapse. The mass 
density, temperature and magnetic field on the progenitor 
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stellar surface can be inferred from the observed range. From 
these quantities, one estimates the dimensionless magnetic 
parameter h, which plays an important role in the MHD 
cvolution of quasi-magnetostatic solutions. By self-similar 
MHD transformation Q and the ideal gas law, we have 



Ck 



l-3g/2 



k B T 



p(47rG)T-iG8(3n - 2)ip~- 1 M<i 



(10) 



where T is the gas temperature, k B is Boltzmann's con- 
stant, and p is the mean molecular (atomic) weight for gas 
particles. For a gas mainly of ionized hydrogen and a stellar 
mass of several solar masses, one can estimate parameter k 
according to expression (j 10|) from 7 and the local density p 
and temperature T. For typical values of p ~ 10~ 5 g cm~ 3 
and T ~ 10 5 — 10 6 K on the stellar surface of a massive pro- 
genitor in the late phase (just before the gravitational core 
collapse) and 7 = 1.3 and n = 0.7, we estimate a range of 
k ~ 10 16 - 10 17 cgs unit. 

With t — > + and/or r — > 00 (i.e. for the regime of 
large x), self-similar MHD solutions follow asymptotic form 

For such analytic asymptotic solutions, the mass density 
simply scales as 



Ak 1/n 
AtvG ' 



-2/n 



(11) 



and is independent of time t. The reason to refer A as the 
mass parameter is now apparent. With asymptotic mass 
density scaling (|ll|l . one can estimate A from the surface 
mass density of a progenitor star. We have the radial fiow 
velocity 



Bk 



V(2n) r l- 



1/n 



(12) 



also independent of t. The radial fiow velocity near the stel- 
lar surface relates to the mass accretion rate or mass loss 
rate by 



M = -iwr pu . (13) 

With expressions (|12[) and (|13[) , one can determine velocity 
parameter B from the stellar mass loss rate for B > or 
mass accretion rate for B < 0. Observationally, mass loss 
rates of Galacti c OB stars fall in the range of ~ 10~ 5 -10~ 7 
Mq yr _1 (e.g. iLamers fc Leithererl 1 19931 ). Mass loss rates 
of Wolf-Rayet stars are much higher than massive main- 
sequence stars, falling in the range of ~ 10~ 4 — 10~ 6 Mp 
yr -1 (e.g. ISinghlll98q ). Very recently, iPuls. Vink fc Naiarrol 
(2008) provide an extensive review on mass loss rates of 
massive stars. 

In summary, with the mass density, temperature, mag- 
netic field and mass loss rate specified at the surface of a 
massive progenitor star, one can determine all parameters of 
asymptotic solution (J5J) and integrate nonlinear MHD ODEs 
([2| and ([3]) inwards to produce an MHD profile for the pro- 
genitor interior and envelope. We still have the freedom to 
choose the rebound shock position after the gravitational 
core collapse. Physically, the speed of the rebound shock de- 
pends on the core collapse, for example the EoS and the neu- 
trino reheating, etc. Here we simply treat it as an adjustable 
parameter to search for downstream solutions within the 
shock front as x — > + . 

With t — > 00 and/or r — * + (i.e. for the regime of 
small x), the final evolution may approach either a quasi- 
magnetostatic solution (Lou & Wang 2007) or the strong- 



field solutions (Wang & Lou 2007, 2008). It is desirable that 
after a long lapse (i.e. t — > 00), the enclosed mass within n 
evolves towards a constant value and form a compact rem- 
nant of nuclear mass density. Asymptotically, the enclosed 
mass takes the form of 

nk 1/n A 3 _ 2/n 



M = 



(3n - 2)G 



(14) 



for quasi-magnetostatic solutions and becomes independent 
of t. This appears consistent with the scenario of form- 
ing a central compa ct object with a strong magnetic field 
jLou fc Wandl2007l ). In a companion paper, we will show 
that the asymptotic enclosed mass for strong-field solutions 
depends on t, and may be invoked to model a continuous 
accretion or outfiow around a nascent neutron star (Hu & 
Lou 2009, in preparation) . As quasi-magnetostatic solutions 
require that h < h c , we arrive at an interesting situation, i.e. 
in order to give birth to a stable neutron star with an ultra- 
intense magnetic field, the massive progenitor star needs to 
be magnetized but not too much. We shall see presently that 
h < h c is generally satisfied for massive main-sequence stars. 
We now introduce the outer initial mass M Q m j and 



the inner ultimate mass M: 



in the same manner as done 



in ILou fc Wand (|2ood I2007I ) and regard them as rough 
estimates for the the initial progenitor mass and the fi- 
nal mass of the remnant compact object, respectively. The 
ratio of these two masses is given by M Q m j /Mj „1+ = 

A*(ro/n) (3 - 2/n) with A* = (A/A )\- 2/n . The ratio of outer 
initial magnetic field at the surface of the progenitor star to 
the inner final magnetic field of the central compact rem- 
nant is < B^ > 1/2 / < B?^ > 1/2 = X* (r / n )^ 2/n K 
The factor A* is insignificant as compared to the radial vari- 
ation of magnetic field, i.e. the very radial dependence of 
r (i-2/«) As n — > 2/3 and the polytropic index 7 approaches 
4/3, this scaling approaches r~ 2 . With ro = 10 12 cm and 
r\ = 10 6 cm, the magnetic field strength can be rapidly en- 
hanced by a factor up to ~ 10 12 . Thus, for a magnetar (i.e. 
neutron star) to have a surface magnetic field strength of 



< B 



i,ult 



> 



1/2, 



10 1S G, we need a magnetic field of ~ 10 3 G 



over the progenitor stellar surface, which is attainable for 
magnetic massive OB stars. 



3.2 Numerical Model Calculations 

The initial time to apply our solutions is estimated by 
the time when the rebound shock crosses n, viz. t\ = 
[n/(k 1/2 x a )} 1/n l|Lou fc Wang||2006l ). Here we suppose that 
roughly from t\ on the collapse and rebound shock inside 
the progenitor star have already evolved into a self-similar 
phase (typically this process takes a few milliseconds) . The 
rebound shock travels outwards throug h the stellar interior 
into the e nvelop e in ~ 10 4 - 10 6 s (e.g. iLou fc Wanell2007l ; 
iLou fc Hu|[2009l ; iHu fc Loul l2008al l . We set t 2 as the time 
when the rebound shock reaches r , roughly around the 
shock breakout. We will also show self-similar MHD solu- 
tions at t m \ = 1 s as an intermediate time between t\ and 
t 2 , and at t — 00. 

As an example of illustration, we choose n = 0.673, q = 
0, 7 = 1.327 for a conventional polytropic gas. From the 
analysis in subsection 3.1, solutions with a smaller n (> 2/3) 
tend to give a larger Mj u jj. and < B 2 U ^ > 1/2 - The choice of 
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Figure 1. The radial profiles of mass density, radial flow velocity, enclosed mass and transverse magnetic field at four different epochs. 
For all panels, the dash-dotted, solid, dashed and dotted curves are at t\, t m \, t^, and t = oo, respectively. Mass density, velocity and 
magnetic field profiles are multiplied by numerical factors marked along the curves for a compact presentation. 



parameters here is a compromise among multiple numerical 
tests. With a typical surface mass density of 2.5 x 10~ 5 g 
cm -3 , a surface temperature 3 x 10 4 K, a mass loss rate 10 -6 
Mq yr _1 , and a surface magnetic field strength 10 3 G for a 
massive progenitor star, we estimate ki = 1.55 x 10 16 cgs 
unit, A = 8.4378, B = 1.27xl0~ 7 and h = 1.52 xlO~ 4 . Such 
parameter h ensures that a stellar core collapse evolves into 
a quasi-magnetostatic manner for small x. Here, inequality 
h < h c is readily satisfied for this set of adopted parameters. 
We set a rebound shock reaching a radius r = 10 cm at 
t = 1 s, which fixes the shock location and travel speed. 
With these parameters, a global self-similar rebound shock 
solution can be constructed and the temporal evolution is 
shown in Figure [JJ at sampled epochs ti = 3.49 x 10 -5 s, 
t m i, t^ = 2.87 x 10 4 s and t = oo. 

The initial progenito r mass is around 5.59 M q, consis- 
tent with observations of ISoderberg et all ([2008) . At epoch 
t\ (dash-dotted curve), the rebound shock has not yet 
emerged. The radial velocities around the surface of the 
proto-neutron star point inwards, corresponding to a core 
collapse process leading to a subsequent rebound shock trav- 
elling outwards inside the progenitor envelope. Meanwhile, 
the outer part of the progenitor envelope still flows slowly 
outwards for a stellar wind. Such kind of self-similar shock 
manifestation with a collapsing inner part and an expand- 
ing outer part is made po ssible from the quas i-magnetostatic 
asymptotic solutions fsee lLou fc W ane 2006) and is another 
form of envelope expansion with core collapse (EECC) pro- 
posed by Lou & Shen (2004). At epoch t m \ (solid curve), 



the outgoing rebound shock emerges and travels within the 
stellar envelope. Figure [JJ clearly shows this discontinuity 
across the outgoing MHD shock front. The rebound shock 
evolves in a self-similar manner with an outgoing speed de- 
creasing with time t for n < 1 (see Lou & Hu 2008 and 
Hu & Lou 2008b for a comparison with numerical simula- 
tions). The immediate downstream side of the shock has 
an outward velocity and the immediate upstream side has 
an inward velocity. The enclosed mass of the downstream 
side decreases rapidly towards the centre while that of the 
upstream side remains nearly unchanged. Across the shock 
front, both mass density and magnetic field strength are en- 
hanced by a factor of 6.96. It can be derived that < Bf >\^ 2 
I < Bl > 1/2 = /n/pa = 2/[( 7 +l)A4?] + ( 7 -l)/( 7 +l) where 
Ali is the upstream Mach number in the comoving shock 
framework of reference. The maximum enhancement across 
the shock is ( 7 + l)/( 7 — 1) = 7.12 for our adopted value of 
polytropic index 7 = 1.327. 

This rebound shock breaks out from the stellar envelope 
m ( 2 ~3x 10 4 s. We see at that moment the flow velocity 
within the spherical volume previously occupied by the pro- 
genitor star becomes very much reduced in the wake of the 
rebound shock, and the gas there gradually approaches the 
quasi-magnetostatic phase of evolution. Coupled with radi- 
ation mechanisms, for instance the thermal bremsstrahlung 
of hot electrons with a temperature T > 10 7 K, and using 
the dynamic profiles shown in Figure [TJ one may compute 
the radiation detected and reproduce the X-ray or 7 -ray 
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Figure 2. The dependence of the ensemble averaged magnetic 

field strength enhancemcnt factor < B? ,, >V 2 / < B 2 . . > 1/2 
° i,ult ' o,im 

and the mass ratio Mj U [ t /M Q ini on thc MHD rebound shock 

radii for n = 0.673, 5 = 0, 7 = 1.327, a surface mass density 

2.5x10 

loss rate 10 -6 M yr" 1 



5 g cm 3 , a surface temperature 3 X 10 4 K, and a mass 



liRht curves observed (e.g. iHu fc Loul2008bl ; ISoderberg et al.l 
120081 ; iMazzali et al.|[2008l V Eventually, flow velocities of the 
entire system tend to zero and the enclosed masses at all 
radii remain unchanged. From Figure[TJ we see that for the 
initial and final stages of the MHD evolution, the mass den- 
sity and magnetic field distributions obey power laws, con- 
sistent with the asymptotic analysis in subsection 3.1. Fi- 
nally, within radius n of the inner compact remnant, the 
enclosed mass is 2.15 Mq with a mean density of 1.02 x 10 15 



g cm 



for a neutron star. As the mean surface magnetic 



field strength is < B? ^ > 1/2 ~ 4.70 x 10 14 G, this neutron 
star should be indeed regarded as a magnetar. 

Numerical explorations indicate that the ultimate mag- 
netic field on the neutron star is proportional to the initial 
magnetic field on the progenitor stellar surface. However, the 
magnetic enhancement factor < -B 2 mt > 1/2 / < B 2 — > 1/2 
and the mass ratio M- x n n/M Q j m indeed depend on model 
parameters and in particular, on shock properties and self- 
similar scaling indices n and q. As long as the density scales 
as r~ 2 ' n , index n must be set to approach the limiting value 
2/3 to ensure a sufficiently massive progenitor star. In Figure 
[2] we plot these two ratios versus the selected rebound shock 
radii at t = 1 s. The two curves suggest that a shock with 
a medium travel speed is associated with a minimum mag- 
netic field enhancement. The magnetic enhancement factor 
appears grossly proportional to the mass ratio. We illustrate 
the relation between these two ratios in FigureO Note that 
the magnetic enhancement factor does not exceed 10 12 . 

We further examine the influence of parameter q for 
a general polytropic magnetofluid shown in Fig. U With a 
larger q, both the mass ratio and magnetic enhancement 
factor become less. The qualitative behaviours of quasi- 
magnetostatic solutions in general polytropic cases are sim- 
ilar to the conventional polytropic cases. 

Our self-similar MHD rebound shock model analysis 
suggests that should there be a continuum of stellar sur- 
face magnetic field strengths over massive progenitor stars, 
there is then a corresponding continuum from normal radio 
pulsars to magnetars in terms of magnetic field strengths 
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Figure 3. The ensemble averaged magnetic field strength en- 

hancement < B? ,, > 1/2 / < B 2 . . > 1/2 versus the mass ra- 
i,ult ' o,im 

tio Afj uit /M Q m j after a long time of rebound shock evolution. 

The star points mark our numerical explorations with different 

shock properties in the case of n = 0.673, q = 0, 7 = 1.327, a 

surface mass density 2.5 X 10 — 5 g cm~ 3 , a surface temperature 

3 X 10 4 K and a mass loss rate 10 6 Mg yr l . Our numerical 

data can be best fitted (solid curve) with a quadratic rclation 

< B lu\t > 1/2 / < B o,ini > 1/2 = ~ 8 - 8 X 10ll ( A/ i,ultA'o,ini) 2 + 
1.8 x 10 12 (Mj lt /M ' ini ) - 7.1 x 10 10 . 
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Figure 4. Thc dcpcndcnce of the ensemble averaged magnetic 

field strength enhancement factor < B ? ^ > 1/2 / < B 2 . . > 1/2 
& i,ult ' 0,1111 

and the mass ratio M^ U j t /M Q ini on scaling index parameter q in 
the case of n = 0.673, a surface mass density 2.5 X 10 — 5 g cm — 3 , 
a surface temperaturc 3 X 10 4 K, a mass loss rate 10 — 6 Mq yr — 1 
and a shock position r s = 7 X 10 8 cm at t = 1 s. 



over the compact stellar surface. The key factors that decide 
whether a remnant neutron star possesses an ultra-intense 
magnetic field include scaling indices n and q, the initial sur- 
face magnetic field of progenitor star, and the strength (or 
speed) of the rebound shock. Ou r analysis al s o resol ves the 
difficulty posed by the survey of iPetit et ahl (|2008bl ). That 
is, a massive progenitor star with a proper range of surface 
magnetic field strengths is most likely but not necessarily 
leads to a magnetar after the gravitational core collapse and 
the emergence of a rebound shock. There are additional con- 
straints to be satisfied. Both Figures [2] and [4] demonstrate 
that the magnetic enhancement factor cannot exceed ~ 10 12 , 
and may be as low as ~ 10 11 even asn-t 2/3. This implies 
that the conditions for producing magnetars are fairly strict. 
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Another major limit arises as the fact that the magnetic en- 
hancement factor correlates to the mass ratio. If this mass 
ratio is too high, the compact remnant mass may exceed the 
Tolman-Oppenheimer-Volkoff (TOV) limit (~ 3 - 3.2 M ; 
Rhoades & Ruffini 1974), and the core object would collapse 
further to form a black hole. The upper bound for the mass 
of neutron stars then places a limit on their surface magnetic 
field strengths. 



4 CONCLUSIONS AND DISCUSSION 

We combine semi-analytic and numerical self-similar MHD 
solutions to model the gravitational core collapse, the re- 
bound shock explosion of a magnetized massive progenitor 
star, and the formation of a central remnant compact mag- 
netar. As natural extension and generalization to recent su- 
pernova rebound shock models of Lou & Wang (2006, 2007), 
we invoked quasi-magnetostatic asymptotic solutions and 
asymptotic solutions far away from the centre for a general 
polytropic magnetofluid. With the magnetic frozen-in condi- 
tion imposed, the surface magnetic field of a nascent neutron 
star can be very much stronger than that of its progenitor 
by a factor of ~ 10 11 — 10 12 during processes of the grav- 
itational core collapse and the self-similar rebound shock 
breakout. Therefore if the progenitor is a magnetic massive 
star with a surface magnetic field strength of ~ 10 3 G, it 
would have a good chance to produce a magnetar at the 
centre of its supernova remnant. Here we propose that mag- 
netars may be produced through powerful supernova explo- 
sions of magnetized massive progenitor stars. Such physical 
origin is also supported by statistical inferences from ob - 
servational surveys (e.g. iFerrario fc Wickramasinghel l2006). 
While the magnetic flux of the collapsing core inside the 
progenitor star may come from either main-sequence stellar 
dynamo processes or 'fossil fields' of molecular clouds, it will 
be dragged and squeezed into the newborn neutron star by 
the conservation of magnetic fluxes. At least for magnetic 
massive stars as magnetar progenitors, the post-supernova 
dynamo processes inside the remnant neutron star may not 
be necessary for producing the ultra-intense surface mag- 
netic field. Magnetic field strengths in the interior of such 
magnetars thus produced are expected to be even stronger 
and are gravitationally buried and confined by the nuclear- 
density matter. 

If the surface magnetic field strength of a massive pro- 
genitor happens to be even higher, possibly up to ~ 3 x 10 4 
G or more, a nascent compact magnetar may then possess 
an ultra-intense surface magnetic field up to ~ 10 15 — 3x 10 16 
G or higher. Such strong magnetic fields can give rise to var- 
ious stellar activities and magnetic reconnection can release 
stored magnetic energy sporadically and violently (e.g. Low 
& Lou 1990). For example, if one approximates the mag- 
netospheres of SGRs, AXPs, and radio pulsars as in gross 
force-free equilibria (i.e. electric currents parallel to mag- 
netic fields), then the magnetic energies retained in such 
magnetospheric systems are higher than the corresponding 
potential field configurations with the same footpoint mag- 
netic field at the stellar surface. A loss of equilibrium most 
likely triggered by magnetic reconnections may then lead to 
outbursts of available magnetic energies. In the solar and 
stellar contexts, such processes correspond to solar/stellar 



flare activities and coronal mass ejections. For compact stel- 
lar object like neutron stars, this type of dramatic magnetic 
energy releases might fuel 'magnetic fireballs' which pro- 
duce short-hard gamma-ray bursts, such as those reported 
giant flares (e.g. iHurlev et alj l2005ħ as well as very recent 
outbursts of SGR J1550-5418 and SGR 1627-41. After gi- 
ant flares, a neutron star may still recover surface magnetic 
fields of order of ~ 10 14 — 10 15 G manifested as either an 
AXP or a SGR. 

For an intense magnetic field buried inside a magne- 
tar, there is yet another possible source of activities. The 
magnetic Lorentz force may break the crust of a neutron 
star into pieces to various levels. Coupled with the mag- 
netar spin, a broken crust can give rise to various activi- 
ties. For examples, chunks of crust may pile up around the 
equatorial bulge and rearrange themselves to generate stel- 
lar seismic activities in a random manner; interior magnetic 
flux ropes may burst into the magnetar magnetosphere ran- 
domly at weak points of a crust; if the crust is more or less 
destroyed by the Lorentz force, then footpoints of magnetic 
field lines can be moved around by convective motions and 
possible differential rotations of a magnetar; the manoeuvre 
of magnetic footpoints over the surface of a magnetar leads 
to variable configurations of the magnetosphere and thus 
produces magnetic activities including analogues of 'flares' 
and 'coronal mass ejections' mentioned earlier. 

Regarding observed magnetar-like X-ray emissions from 
a rotation-powered pul sar PSR J1846-0258 inside supernova 
remn ant Kes75 (e.g. iGavriil et al. 2008; lArchibald et alj 
2008), our model can accommodate this pulsar resulting 
from a magnetized massive progenitor yet with a lower sur- 
face magnetic field strength. Such magnetar-like activities 
are physically associated more with strong magnetospheric 
field strengths. With the current observational evidence, it 
appears not necessary to postulate that magnetars evolve 
from fast spinning radio pulsars. 

In our 'ball of thread' scenario for a random magnetic 
field within a massive progenitor star, the large-scale mean 
of such a random magnetic field is idealized as dominantly 
transverse with a fairly weak radial component. By the 
approximation of quasi-spherical symmetry, low-amplitude 
small-scale deviations, oscillations or fluctuations are ran- 
domly distributed about the mean flow profile and are ex- 
pected to co-evolve with large-scale MHD profiles (e.g. Lou 
& Bai 2009 in preparation; Cao & Lou 2009 in preparation) . 
During the processes of the gravitational core collapse and 
the MHD rebound shock breakout, a remnant neutron star 
forms with high nuclear density and ultra-intense magnetic 
field, while a major portion of the interior and envelope of 
the massive progenitor star is driven out into the interstellar 
space by the rebound shock with entrained strong magnetic 
field. Our semi-analytic model describes a large-scale self- 
similar MHD evolution for the supernova explosion of a mag- 
netized massive progenitor. Along with this large-scale evo- 
lution, the central magnetar and the thrown-out stellar ma- 
terials may certainly follow their own courses of adjustment 
or rearrangement (e.g. Lou 1994). For example, considerable 
radial components of magnetic field should emerge in a ran- 
dom turbulent manner. In fact, magnetars are expected to 
possess magnetospheres with various possible configurations 
(see Low & Lou 1990 for constructing force-free stellar mag- 
netic field configurations) . By numcrical simulations with 
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'fossil fields'. iBraithwaite fc SpruitJ (|2004I ) illustrated exam- 
ples of processes for magnetic field rearrangement to occur 
within a few Alfven timescales. The emergent stable struc- 
ture of magnetic field for magnetic Ap stars appears always 
of 'offset dipole' type (with complex and twisted magnetic 
field configurations inside), consistent with observations. By 
this analogy, we presume that such magnetic field rearrange- 
ment processes could take place very rapidly during the for- 
mation of intensely magnetized neutron stars, whose Alfven 
timescale is of the order of ~ 0.1 s. Eventually, a magnetar 
can possess a variety of magnetic field configurations (e.g. 
Low & Lou 1990). 

Magnetars observed so far appear to be slow rotators, 
while massive stars are in general rapid rotators with typi- 
cal equatorial speeds of ~ 200 km s" 1 (e.g. Fukuda 1982). 
Therefore, significant angular momentum transfer may have 
taken place i n the stella r evolution of magnetic massive 
stars. ISpruitl l| 19991 . 12002? ) has shown that magnetic fields 
can be created in stably stratified layers inside a differen- 
tially rotating star. Heger et al. (2005) gave detailed rotating 
stellar evolution calculations for stars in the mass range of 
~ 12 — 35 Mq incorporating the dynamo-powered magnetic 
field. In general, it is found that magnetic breaking decreases 
the final spin rate of the collapsing iron core by a factor of 
~ 30 — 50 when compared with the nonmagnetic case. The 
'fossil' (or primordial) magnetic fields may have similar dy- 
namic effects regarding the re-distribution of angular mo- 
mentum inside a massive star. In particular, for magnetar 
formation, high magnetic fields may lead to stronger core- 
envelope coupling during the hydrogen and helium burning 
phase of the SN progenitor, and the collapsing iron core 
and the compact supernova remnant is expected to be even 
slower rotators. This is in accordance with the population 
of slowly rotating magnetars observed. 

In our model at this stage of development, the stel- 
lar rotation is not included to simplify the mathematical 
treatment. Conceptually, it could be possible to design an 
axisymmetric MHD problem to accommodate stellar differ- 
ential rotation in order to explore the re-distribution of an- 
gular momentum during the processes of gravitational core 
collapse, MHD rebound shock as well as collimated MHD 
outflows from polar regions with shocks (e.g. shocked MHD 
jets). The overall magnetic field configuration could be pre- 
dominantly toroidal but a relatively weak radial magnetic 
field component is necessary to exert an effective magnetic 
torque to break or slow down the stellar core rotation. As 
the core materials rapidly collapse towards the centre un- 
der the gravity, the mechanical angular momentum is trans- 
ferred outwards in an outgoing envelope with shock. Along 
the rotation axis, collimated outflows or jets may emerge 
to breakthrough the polar stellar envelope and part of the 
mechanical angular momentum is carried outwards by ro- 
tating polar collimated outflows. For a semi-analytic self- 
similar approach to this time-dependent problem, one might 
be able to perform a self-similar transformation combining 
time t with two spatial coordinates, say r and 8. It might 
be possible to derive asymptotic solutions in the regime of 
slow rotators for this two-dimensional magnetar formation 
problem. One also expects the existence of several MHD sin- 
gular surfaces in deriving flow solutions. Physically, such a 
scheme if tractable semi-analytically and/or numerically can 
be applied to a wide range of gravitational collapses of rotat- 



ing systems, including magnetars, pulsars, magnetic whitc 
dwarfs, protostars, planets and so forth. 

In addition to the quasi-magnetostatic asymptotic so- 
lutions adopted and exemplified in this paper, it may be 
possible that a magnetized massive progenitor star evolves 
towards the strong-field asymptotic sol utions ultimately 
llYu et al. I l200d : IWang fc Loull2007l , l2008f l. involving a ma- 
terial fall-back process towards the central remnant neutron 
star. Under this situation, the enclosed mass within a certain 
radius keeps increasing until the mass of the neutron star ex- 
ceeds the TOV limit. Such an MHD fall-back process offers 
a possible means to form stellar mass black holes as compact 
remnants in supernova and hypernova explosions. We em- 
phasize that such a mechanism requires a strong magnetic 
field inside a progenitor and the magnetic force becomes 
dominant during the fall-back process. 
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jWang fe Loui r2008'). Once we have (x^, Q2, «2) on the im- 
mediate downstream side of a shock front, we can obtain 
(xi, Qi, v\) explicitly for the immediat e upstream side us- 
ing MHD shock relations (JXT) - (O l|Wang fe Loul |2008| ; 
iLou fc Hul l2009h or vice versa. In the case of q = 2/3, 
there are only two independent relations among equations 
(| All) — (| A3() and we could choose parameter A > arbitrar- 
ily. Hence we can set fci = or A = 1 in this situation. 
This treatment will not alter the relations of the result- 
ing dimensional physical variables. In general, the outgoing 
travel speed of a rebound shock varies with time t for n 7^ 1: 
shock accelerates for n > 1, shock speed remains constant 
for n — 1, and shock decelerates for n < 1. 



APPENDIX A: MHD SHOCK CONDITIONS 

In the shock comoving framework of reference, the MHD 
shock jump conditions in terms of the reduced self-similar 



